
*** Do File for: Authoritarianism and Support for Social Security PRQ ***

*Regressions
clear 
use "C:\Users\dmac9\Dropbox\Accepted Authoritarianism and Social Security (PRQ)\CANES Data.dta"

*CANES, full sample*
probit more_ss Auth_Index01 female white married own_home Age62 Retired Disabled income01 College Confed pid7_01 ideo7_01 i.year [pweight=weight]
*outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') replace
outreg2 using tab.rtf, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') replace
margins, at(Auth_Index01=(0(.25)1)) 
marginsplot

*Additional Programs*
clear 
use "C:\Users\dmac9\Dropbox\Accepted Authoritarianism and Social Security (PRQ)\ANES 2016 Data.dta"

probit more_ss auth_index01 female white age age2 married own_home Income_cat01 college Confed Pid7_01 Ideo7_01 online_mode [pweight=weight] 
outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') replace
margins, dydx(auth_index01) saving(m1, replace)

probit more_childcare auth_index01 female white age age2 married own_home Income_cat01 college Confed Pid7_01 Ideo7_01 online_mode [pweight=weight]
outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') append
margins, dydx(auth_index01) saving(m2, replace)

probit more_poor auth_index01 female white age age2 married own_home Income_cat01 college Confed Pid7_01 Ideo7_01 online_mode [pweight=weight]
outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') append
margins, dydx(auth_index01) saving(m3, replace)

probit more_schools auth_index01 female white age age2 married own_home Income_cat01 college Confed Pid7_01 Ideo7_01 online_mode [pweight=weight]
outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') append
margins, dydx(auth_index01) saving(m4, replace)

probit more_welfare auth_index01 female white age age2 married own_home Income_cat01 college Confed Pid7_01 Ideo7_01 online_mode [pweight=weight]
outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') append
margins, dydx(auth_index01) saving(m5, replace)

probit more_health_spending auth_index01 female white age age2 married own_home Income_cat01 college Confed Pid7_01 Ideo7_01 online_mode [pweight=weight]
outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') append
margins, dydx(auth_index01) saving(m6, replace)

probit tax_millionaires auth_index01 female white age age2 married own_home Income_cat01 college Confed Pid7_01 Ideo7_01 online_mode [pweight=weight]
outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') append
margins, dydx(auth_index01) saving(m7, replace)

probit raise_minwage auth_index01 female white age age2 married own_home Income_cat01 college Confed Pid7_01 Ideo7_01 online_mode [pweight=weight]
outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') append
margins, dydx(auth_index01) saving(m8, replace)

probit more_crime auth_index01 female white age age2 married own_home Income_cat01 college Confed Pid7_01 Ideo7_01 online_mode [pweight=weight]
outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') append
margins, dydx(auth_index01) saving(m9, replace)

probit more_defense_spending auth_index01 female white age age2 married own_home Income_cat01 college Confed Pid7_01 Ideo7_01 online_mode [pweight=weight]
outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') append
margins, dydx(auth_index01) saving(m10, replace)

combomarginsplot m1 m2 m3 m4 m5 m6 m7 m8 m9 m10, ///
labels("Social Security" "Child Care" "Aid to Poor" "Public Schools" "Welfare Programs" "Health Care" "Tax Millionaires" "Minimum Wage" "Dealing with Crime" "Defense Spending") horizontal


*Bivariate vs. Simpler model vs. More Complex model ANES*
clear 
use "C:\Users\dmac9\Dropbox\Accepted Authoritarianism and Social Security (PRQ)\CANES Data.dta"

probit more_ss Auth_Index01 [pweight=weight] 
outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') replace
margins, dydx(Auth_Index01) saving(m1, replace)

probit more_ss Auth_Index01 female white married own_home Age62 Retired Disabled income01 College Confed pid7_01 ideo7_01 i.year [pweight=weight] 
outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') append
margins, dydx(Auth_Index01) saving(m2, replace)

probit more_ss Auth_Index01 female white married own_home Age62 Retired Disabled income01 College Confed pid7_01 ideo7_01 ///
egal_index01 lg_index01 mt_index01 ft_poor01 whites_hardworking01 blacks_hardworking01 hispanics_hardworking01 asians_hardworking01 i.year [pweight=weight]
outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') append
margins, dydx(Auth_Index01) saving(m3, replace)

combomarginsplot m1 m2 m3, labels("No Controls" "Main Model" "More Controls") 


*GSS panels*
clear 
use "C:\Users\dmac9\Dropbox\Accepted Authoritarianism and Social Security (PRQ)\GSS Panel 08-12 Data.dta"

reg Obey_2012_01 Obey_2008_01 More_SS_2008, r
*outreg2 using tab.tex, dec(3) replace
outreg2 using tab.rtf, dec(3) replace
reg More_SS_2012 Obey_2008_01 More_SS_2008, r
*outreg2 using tab.tex, dec(3) append
outreg2 using tab.rtf, dec(3) append

clear 
use "C:\Users\dmac9\Dropbox\Accepted Authoritarianism and Social Security (PRQ)\GSS Panel 10-14 Data.dta"

reg Obey_2014_01 Obey_2010_01 More_SS_2010, r
*outreg2 using tab.tex, dec(3) append
outreg2 using tab.rtf, dec(3) append
reg More_SS_2014 Obey_2010_01 More_SS_2010, r
*outreg2 using tab.tex, dec(3)  append
outreg2 using tab.rtf, dec(3) append

*Heterogeneity in the Relationship*
clear 
use "C:\Users\dmac9\Dropbox\Accepted Authoritarianism and Social Security (PRQ)\CANES Data.dta"

*By Race/Ethnicity*
probit more_ss Auth_Index01 female white married own_home Age62 Retired Disabled income01 College Confed pid7_01 ideo7_01 i.year [pweight=weight]
outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') replace
margins, dydx(Auth_Index01) saving(m1, replace)

probit more_ss Auth_Index01 female married own_home Age62 Retired Disabled income01 College Confed pid7_01 ideo7_01 i.year [pweight=weight] if white==1
outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') append
margins, dydx(Auth_Index01) saving(m2, replace)

probit more_ss Auth_Index01 female married own_home Age62 Retired Disabled income01 College Confed pid7_01 ideo7_01 i.year [pweight=weight] if white==0
outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') append
margins, dydx(Auth_Index01) saving(m3, replace)

combomarginsplot m1 m2 m3, labels("Full Sample" "Whites" "Non-Whites") 

*Political Information*
clear 
use "C:\Users\dmac9\Dropbox\Accepted Authoritarianism and Social Security (PRQ)\CANES Data.dta"

probit more_ss Auth_Index01 female white married own_home Age62 Retired Disabled income01 College Confed pid7_01 ideo7_01 i.year [pweight=weight] if Int_pk_pre3==1
outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') replace
margins, dydx(Auth_Index01) saving(m1, replace)

probit more_ss Auth_Index01 female white married own_home Age62 Retired Disabled income01 College Confed pid7_01 ideo7_01 i.year [pweight=weight] if Int_pk_pre3==2
outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') append
margins, dydx(Auth_Index01) saving(m2, replace)

probit more_ss Auth_Index01 female white married own_home Age62 Retired Disabled income01 College Confed pid7_01 ideo7_01 i.year [pweight=weight] if Int_pk_pre3==3 
outreg2 using tab.tex, dec(3) addstat(Pseudo R$^2$, `e(r2_p)') append
margins, dydx(Auth_Index01) saving(m3, replace)

combomarginsplot m1 m2 m3, labels("Low Info" "Avg Info" "High Info") 


